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Abstract 

We construct the phase diagram of the quark-antiquark and diquark condensates at finite tem- 
perature and density in the 2 + 1 dimensional (3D) two flavor massless Gross Neveu (GN) model 
with the 4-component quarks. In contrast to the case of the 2-component quarks, there appears 
the coexisting phase of the quark-antiquark and diquark condensates. This is the crucial difference 
between the 2-component and 4-component quark cases in the 3D GN model. The coexisting phase 
is also seen in the 4D Nambu Jona-Lasinio model. Then we see that the 3D GN model with the 
4-component quarks bears closer resemblance to the 4D Nambu Jona-Lasinio model. 
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I. INTRODUCTION 



In recent years, a lot of works has been devoted to the study of the phase structure of 
Quantum Chromodynamics (QCD). QCD is an asymptotically free theory and the interac- 
tions between quarks and gluons become weak at high energy 1] . Then, at high temperature 
and/or density, the quarks and gluons constitute rather weakly interacting system, which is 
called the quark-gluon plasma. On the other hand, at low temperature and density, quarks 
and gluons are confined into hadrons and can not be observed as free particles. Further- 
more, the existence of the color superconducting phases at low temperature and moderate 
baryon density has been widely accepted. The color superconductivity is the state where 
the quark-quark (diquark) Cooper pairs [2j are induced by the attractive interaction in the 
color antitriplet channel[3|. 

One of the simple theories to study the above mentioned subject is the Nambu Jona- 
Lasinio (NJL) model, which is a low energy effective field theory of QCD [4]. The NJL 
model successfully describes the QCD phase structure, and a variety of works has been 
devoted to the study on the basis of the NJL model. (For nice reviews, see, e.g. {5, 6].) In 
particular, through analysis on the competition between quark-antiquark (qq) and diquark 
(qq) condensates in the NJL model, it has been found that there appears the coexisting 
phase of the qq and qq condensed phases (see js]). 

The study of the phase structure of the NJL model in lower dimensions (D) is also an 
interesting issue since the models usually become simpler in lower dimensions 7|] . Indeed, 
the Gross Neveu (GN) model 8|, which is the counterpart of the NJL model in 2D and 
3D, becomes renormalizable. Recently, the phase structure of the qq and qq condensates 
was studied within the 3D GN modeip, Q. There the quarks are assigned to the lowest 
nontrivial (2-dimensional) representation of the 0(2, 1) group which we refer to as the 2- 
component (2c) quarks. The resultant phase diagrams show that there does not appear the 
region where the qq and qq condensates coexist, which is in sharp contrast to the 4D NJL 
case. This difference may stem from the difference in form of the qq condensate term in 
the Lagrangian density between both cases. In the 2c quark case in 3D, the qq condensate 
term in the Lagrangian does not include the 7 s matrix, which leads to the above-mentioned 



difference. However, as discussed in 



111 ], the 7 5 enters in the case of the 4c quarks in 3D. In 



this sense, the 3D GN model with the 4c quarks is expected to bear closer resemblance to 
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the 4D NJL model. Actually, in vacuum (zero temperature and chemical potential) theory, 
the coexisting phase of the qq and qq condensates has been found in the 3D GN model with 



the 4c quarks 



121. 



In this paper, we study the qq and qq condensates at finite temperature and density in 
the 3D 2 flavor massless GN model with the 4c quarks. We construct the phase diagrams, 
and discuss the similarities and differences among the models, the 3D GN model with the 
4c quarks, the 3D GN model with the 2c quarks and the 4D NJL model. 

The paper is organized as follows: In Sec. II we introduce the Lagrangian density in the 
3D GN model with the 4c quarks and employ the mean-field approximation. In Sec. Ill the 
derivation of the thermodynamic potential is given. Then, in Sec. IV, we show the numerical 
results for the qq and qq condensates at zero and finite temperature. In Sec.V, the phase 
diagrams of the qq and qq condensates are obtained. Sec. VI is devoted to concluding remarks. 



II. 3D GROSS NEVEU MODEL 



Following the same reasoning as described in 6], we employ the Lagrangian density: 

C=qi$q + G s (qq) 2 



+ G D (qiT2\2lbq C ){q C iT2^2lbq)- 



(1) 



The notations used in Eq.flT]) are the same as in the 3D 2 flavor massless GN model with 
the 2c quarks given in [9(, except that the 7 matrices are 4x4 matrices here. For the 7 



matrices, we use the same form as in 



7 



7 



73 W=[ iTx 

-r 3 J ' \o -in 

iT2 V- J 1 

-it 2 j ^-10 

The charge conjugated fields are given by q c = Cq T , q c = q T C, where C = 7 2 . 
Employing the mean-field approximation, we have 

1 



(2) 



C =q(i0 - a)q + -A*(q c tr 2 X 2l5 q) 



+ ^ A (^ T 2 A 2759 C ) 



a 
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D 



(3) 
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Here a and A are the order parameters for the qq and qq condensates, which are defined by 

a = -2G s (qq) , A = 2G D (q C ir 2 \ 2l5 q) . (4) 

III. THE THERMODYNAMIC POTENTIAL 

The partition function of the ground canonical ensemble is calculated by using the stan- 
dard method, 

Z = N' J [dq][dq] ex p{^r Jd 2 *(C + m>?)}, (5) 

where /i is the quark chemical potential and j3 = 1/T is the inverse temperature. Introducing 
the Nambu-Gorkov basis [3] 

*= ( \ I and * = (q f) , (6) 



q c 



we have 



Z =N' exp | - jf dr J d 2 x ^ 



a 2 |A| 2 
+ 



4Gc 4G 



D 



hd^exp^-^iPG- 1 )^. (7) 



The matrix G 1 is defined by 



ir 2 A 2 7 5 A* (^-o--/i7 )l/l c 

where 1 / and l c are the unit matrix in flavor and color spaces, respectively. One can compute 
the thermodynamic potential Vt = — \nZ / (3V by following the same procedure as in [6J, 

n(a,\A\) = n (a,\A\)+n T (a,\A\), (9) 
a 2 |A| 2 r A2 n 



no(<r,|A|) 



4G 5 4G Z 

(2^)2 _ 



ln(l + e- /JB± )+21n(l + e-^; 



(11) 



+ 4^- 4 /(^^ + ^ + + ^' (10) 
Sl T (a,\A\) = -4Tj2 [■ 

where V is the volume of the system and 

E^ 2 = {E± jj,) 2 + |A| 2 , E ± = E±/2, 

E = y/p 2 + a 2 , p 2 =p 2 +p 2 . (12) 
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Oo is T independent contribution and VLt is the T dependent part. VLq is ultraviolet divergent 
and we carry out the renormalization through introducing the counter Lagrangian density^, 
which is of the form Cc = —Zs<y 2 /2 — Zd\A\ 2 . In the present case, the renormalization 



factors Zs and are given by 



Z s 
Z D 



12 . 3 
— A a 

7T 2 
4 A 1 

-A a. 

7T 2 



(13) 
(14) 



Here A is the 3D momentum cut-off and a is an arbitrary renormalization scale. 
Then, by introducing the following parameters, 



A 



7T / 1 

3 \4G^ 

7T / 1 

"2 V4G^ 



3 

-a 
4 

1 

r 



we obtain the renormalized thermodynamic potential: 
Q r (a, |A|) = Q 0r (cr, |A|) + fi T ((x, |A|) , 



(15) 
(16) 

(17) 



^OrCc, |A| 



+ 



-E 

3tt 4^ 



3(T 

a 

7T 

(2a 2 + 2|A| 2 - /i 2 ± fjur)y/o* + | A| 2 + /i 2 ± 2fw 



^|A| 2 + -a 3 

7T 37T 



=F 3/i| A| 2 ln{a ± // + ^a 2 + \ A| 2 + ± 2^} 



(18) 



ctq (Aq) is the gg (gg) condensate at T = and /x = in the model when the gg (gg) 



condensate is absent. It shou 
parameters Oq and A . As in 



d be noted that the thermodynamic potential has two free 
^1, we take a to be the scale of the theory (cr > 0). Then, 
after fixing <To, there remains one free parameter A and we study the various values for the 
ratio A /o"o- 



IV. QUARK- ANTIQU ARK AND DIQUARK CONDENSATES 

In this section, we show the numerical result of the gg and gg condensates through 
analyzing the thermodynamic potential, Eq. (fT7|) . The realized condensates are obtained by 
finding the minimum of the thermodynamic potential with respect to a and A. 

Fig. [1] displays the results of the two condensates against /j at T — 0. From the panel 
(a) (A /cr = —6), we see that the gg condensate exists for small /i and it disappears 

5 




FIG. 1: o (circles) and A (triangles) as a function of chemical potential \i at T = 0. 

when \x becomes large. This clearly shows the phenomena of the phase transition, and the 
transition chemical potential is /i = l.Ooo- There does not occur the qq condensate for whole 
/i. However, for A /o"o = —1, the qq condensate appears at /i — l.Ooo, the value where the 
qq condensate disappears. Similar results are obtained in the cases A /o"o = 0, 1/2. At some 
chemical potential, the qq condensate falls and the qq condensate arises. The transition 
chemical potentials are /x = 0.9(T an d 0.8(T for A /cr = and 1/2, respectively. Thus the 
transition chemical potential becomes smaller with increasing the ratio A /o"o! an d the value 
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FIG. 2: The two gaps a (circles) and A (triangles) against T. 

of the qq condensate A enlarges. For A /o"o = 1.1, the qq condensate at fi = has the close 
value with the qq condensate, and it eventually exceeds the qq condensate for A /a = 3/2. 
Note that the qq condensate does not exist for Ao/Vo = 3/2 and there appears only the qq 
condensate. Thus we see that the behaviors of the qq and qq condensates are sensitive to 
the value of the ratio Ao/o~o- 

The results of the condensates at finite temperature are shown in Fig. [3 In the panel (a) 
(A /cr = 0, fi = 0), we see that the qq condensate occurs at T = and it decreases when 
T becomes larger. The transition temperature in this case is T = 0.65o"o. The panel (b) 
(A /o"o = 0, fi = l.Ooo) shows that the qq condensate is 0.65o"o at T = and it disappears 
at T = 0.35a . Thus we see that the qq and qq condensates decreases as the temperature 
increases. This is also the case for Ao/o"o = 1.1, which is shown in the panels (c) and (d). 
We have analyzed the cases for other values of Ao/cxo, and found that the behavior of the 
two condensates do not change qualitatively. 

It should be noted that, in Fig. [T] (T = case), the qq condensate disappears discon- 
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FIG. 3: The phase diagrams. 

tinuously. This is the signal of the first order phase transition. On the other hand, the qq 
condensate disappears continuously in Fig. 12(a) (/i = case), which indicates the second 
order phase transition. This means that the critical point from the first order phase tran- 
sition to the second order one appears at some point in the phase diagram, which we will 
discuss in more detail in the next section. 



V. THE PHASE DIAGRAM 

On the basis of the results obtained in the previous section, we construct the phase 
diagram. Fig. [3] displays the phase structure of the qq and qq condensates for the cases 
AoAxq = -1, 0, 1/2, 1.1. 

From the panel (a), we see that the qq condensate phase is realized at low T and /x, and 
the qq condensate phase occurs at low T and high /x. At high T and /i, there does not exist 
condensate, namely the system is in the normal phase. The similar structure is seen in the 
panel (b): qq phase at low T and /i, qq phase at low T and high [i. However, in the case of 
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Ao/tfo = 0, we see that the coexisting phase of the qq and qq condensates appears at low T 
and intermediate \x. We refer this phase to as the "Double Broken" (DB) phase. Thus, the 
system is characterized by the following phases: 

Q : qq condensate phase [a ^ 0, A = 0) 
D : qq condensate phase (a = 0, A ^ 0) 
DB : Double broken phase (o-^O.A^O) 
N : Normal phase (a = 0, A = 0) 

All of these phases appear in the case of Ao/o~o = 1/2 (panel (c)). However the Q phase 
does not appear for A /o" = 1.1 and only D phase and DB phase appear. Thus, as the ratio 
Ao/o"o increases, the Q condensate phase shrinks toward \i axis and the region of D phase 
and DB phase become larger. 

The points (/i c , T c ) in the panels (a), (b) and (c) indicate the critical points from the 
first order phase transition to the second order one. The phase transition below the critical 
temperature T c is of the first order and above T c is of the second order. We notice that the 
value of T c increases when A /cxo becomes larger, while \x c decreases with increasing A /a . 
Then the critical point moves upward along the transition line, and the region of the first 
order phase transition expands. The other phase transitions, namely the transitions Q — > 
DB, Q — ► N and D — > N, are of the second order. 

VI. CONCLUDING REMARKS 

Through studying the qq and qq condensates at finite temperature and chemical potential, 
we have obtained the phase diagrams in the 3D GN model with the 4c quarks. 

We have shown that the qq condensate (Q) phase is realized at low T and //, the double 
broken (DB) phase at low T and intermediate /i, and the qq condensate (D) phase at low 
T and high /i (see Fig. MJo) and (c)). This feature bears resemblance to the case of the 
4D NJL model. It is difficult to make a direct comparison between the present model and 
the 4D NJL model, because the free parameter in the NJL model is the "direct ratio" of 
the coupling constants Gd/Gs, while the free parameter of the present model is the ratio 
A /o" which is not G^/Gs- However the parameter A is related to Go through Eq. ffToD . 
and Aq becomes larger as Go increases. Then the ratio Aq/cto increases as the qq coupling 
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constant Gd increases. This means that, as Gd/Gs increases, the region of the Q phase 
shrinks, and the regions of the D phase and the DB phase expand. Then, the behavior of 
the phase diagrams shows the close similarity to the case of the 4D NJL model. 

We are now in the position to make the comparison with the 3D GN model with the 2c 
quarks. Comparing with the phase diagrams in the 2c quark case obtained in |9|], we see 
rather similar structure: The Q phase at low T and fi, and the D phase at low T and high 
fi. However, there is one crucial difference. There does not appear the DB phase in the 2c 
quark case, while in the present 4c quark case, the DB phase does appear. The latter fact 
is in accord with the expectation mentioned in the introduction. 

The critical points between the first-order and second-order transition, are located on the 
Q — ► D phase transition line in Fig. E^a) and on the DB — ► D transition line in Fig. 
and (c). In the model without qq condensate, the qq phase transition at zero temperature is 
of the first order and the transition at finite temperature is of the second order. This means 
that T c is negligibly small in the model without the qq condensate. While in the present 
model, T c is finite, which is the reflection of the existence of the qq condensate. It is also 
worth comparing the present model with the model with the 2c quarks, where the coexisting 
phase is absent. The phase transition Q —>■ D is always of the first order and the critical 
points are seen on the line Q -> nQ. This may come from the fact that the existence of the 
qq condensate expels the qq condensate. On the other hand, in the present 4c quark case, 
the qq and qq condensates can coexist and the Q — > D and DB — > D transitions can be of 
the second order. 

Finally, it is worth reemphasizing that the significant qualitative difference is that there 
exists the qq and qq coexisting phase in the present model, which is not seen in the 2c quark 
case. Then, when compared to the case of the 3D GN model with the 2c quarks, the phase 
structures of the present 4c quark case bear closer resemblance to the 4D NJL model. 
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